Levitation and percolation in quantum Hall systems with correlated disorder 
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We investigate tlie integer quantum Hall system in a two dimensional lattice model with spatially 
correlated disorder by using the efficient method to calculate the Chern number proposed by Fukui 
et al. Distribution of charge density indicates that the extended states at the center of each Landau 
band have percolating current paths, which are topologically equivalent to the edge states that exist 
in a system with boundaries. As increasing the strength of disorder, floating feature is observed in 
an averaged Hall conductance as a function of filling factor. Its relation to the observed experiments 
is also discussed. 

PACS numbers: 73.43.-f, 71.30.+h, 71.55.Jv 



Two dimensional electron systems with disorder in a 
perpendicular magnetic field have attracted much atten- 
tion both experimentally and theoretically. According to 
the scaling theorj*i, all single-electron states are local- 
ized at zero temperatures (Anderson localization). How- 
ever, in the strong magnetic field, all states form Landau 
levels and there exist delocalized states at the center of 
each Landau level, which are extended and form per- 
colating paths in a context of network model^i^. The 
transition from the Hall liquid to the Anderson insula- 
tor in the limit of S — > was discussed by Khmelnit- 
skii and Laughlin in terms of the levitation scenario, i.e., 
floating theory^i^, where Landau bands merge together 
as decreasing the magnetic field and the extended states 
float up towards higher energy through the Fermi level. 
This scenario gives a reasonable explanation for many ex- 
periments and the global phase diagram was constructed 
based on it^. However, there still remain many contro- 
versial arguments about the microscopic picture of the 
floating theory^. 

In the quantum Hall system, the topological structure 
of the ground state wave function plays an essential role 
for quantization of the Hall conductance. Each quan- 
tum Hall state can be assigned a topological quantum 
number as the Chern number^, which is the Hall con- 
ductance in unit of e'^ /h. Although there are several 
work a^i^°'^^i^^'^'^ to calculate the Chern numbers and suc- 
ceeded to explain experimentally observed direct transi- 
tion, the random potential, which is spatially uncorre- 
lated, is not compatible with the continuum limi t^^'^^ . 
Koschny et al. proposed that unless the random poten- 
tial has long range correlation, the floating theory in a 
lattice system cannot be consistent with that in a con- 
tinuous system- . By investigating the localization length 
using a recursive Green-function methodic, they found 
the negative-Chern number cannot be created if the cor- 
relation length of randomness is long enough. However, 
since the main part of their study is to calculate Green- 
function, topological nature of extend states is unclear. 

To study topological nature of one-electron states, we 
numerically calculate both the spatial distribution of 
wave functions and the Chern number in this paper. 



Here, we make a spatially correlated random potential in 
order to explore features of a continuous system in a lat- 
tice system. In previous studies on the floating scenario 
with the correlated disorder-i^ii^iiiiiSiii^, the connection 
between the spatial probability density and the Chern 
number has not been focused on. To calculate the Chern 
number, an efflcient numerical method recently proposed 
by Fukui et al.^'^ is employed. 

A tight-binding Hamiltonian on a square lattice with 
nearest neighbor hopping is given by 
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The magnetic flux per plaquette (j) is given by 
J2piaquette ^r,r' = 2TT(f>, whcrc the Summation runs over 
four links around a plaquette. We choose the Landau 
gauge, that is, 9r,r+x = 0, 9r,r+y — 2TT(j)rx. We set the 
units of energy to t — 1, and mostly set the magnetic 
field per one plaquette = 1/16. A random potential, 
Wr, is spatially Gaussian correlated as 



(0.2) 



with mean value (wr) = 0, where ( ) denotes the ensem- 
ble average over different randomness realizations. The 
distribution of the random potential is characterized by 
the correlation length A and the strength W . Although 
the A dependence is importan t ^^i"'^^ , we will focus on the 
dependence on W in the present paper and A is fixed to 
a large value enough to consider a continuous system in 
a lattice system. 

In order to calculate Hall conductance in a system 
without the translational symmetry, we impose a gener- 
alized boundary condition on every randomness realiza- 
tion, by requiring T^(i^)|V^) = exp(i6'^)|V') [[i = x,y), 
where Tfj_{L^) is the magnetic translation operator for a 
displacement and is the boundary-condition phase. 
The system size is L ~ x Ly. The Hall conductance 
axy can be calculated by 6 dependence of the wave func- 
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where is a vector defined as 6 — 



'XI "y 



). As 



for fixed particle density we define a ground state multi- 
plet xjj{0), which is a L x q matrix written as tl}{9) — 
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{\MO)):-;\Mm with H{e)\^,{e)) = e,{e)\i,,)2i. 

Here, g is a number of one-particle states below the Fermi 
energy. Then we define a non-Abelian Berry connection 
^ as ^ = ip'^dif). Finally, the associated Chern number 
is given by 

Cv, = / Tr dA, (0.3) 

where S" is a two dimensional torus parameterized by 
the boundary-condition phases 9^ and 9y. Recently, 
Fukui et alM developed an efficient method to calculate 
a Chern number on a discrete momentum space by ap- 
plying the geometrical formulation of topological charges 
in the lattice gauge theory. As an advantage of their 
method, one can obtain the Chern number accurately 
even if one uses coarsely discretized mesh points in the 
momentum space. 

When randomness is introduced to a pure system, the 
Landau bands will broaden. At the tails of each band, 
all one-electron states become localized due to the An- 
derson localization. On the other hand, there still exist 
extended states at the center of each Landau band. Fig- 
ure [T] shows the dependence of Hall conductance axy and 
density of states (DOS) on the Fermi energy Ef for the 
given randomness realization that is plotted in Fig. [D^a) 
as a contour map. The strength of the randomness W is 
set to 0.01. Nevertheless the energy gap induced by the 
randomness is small, numerical calculations of the Hall 
conductance are stable. As shown in Fig.[Tl the Hall con- 
ductance jumps whenever the Fermi energy goes across 
the center of each Landau band. The extended states 
in the center of each Landau band can be classified into 
two classes; one possess non-zero Chern number, which 
is called a critical state, and the other does not. There 
is only one critical state in the middle of the lowest band 
as well as the second lowest band. In general, there can 
exist several critical states at the center of Landau band. 

To study topological structure of a wave function, 
the contour maps of the probability densities for several 
states are shown in Fig.[2ljb), (c) and (d). Note that we 
adopted the periodic boundary condition. As shown in 
Fig. IHb), the probability density of the localized state 
in the lowest band tail, is concentrated around the valley 
of the random potential. On the other hand, the critical 
state as seen in Fig. [2ljc) or (d) forms a percolating path 
along the y-direction, which is associated with the per- 
colation of equipotential lines2.i^. As shown in Fig.[2fd), 
two percolating paths and one nodal line between them 
are observed for the second Landau bands. This distri- 
bution reminds us that the wave functions in the Landau 
gauge can be described by the Hermite function iJ„(r), 
which has n — 1 nodal lines. Since other extended states 
which do not have nonzero Chern number show the sim- 
ilar spatial distribution, the percolating path is thought 
to be a characteristic feature of extended states. 

There is a relationship between the topological charac- 
ter and the spatial distribution of wave function. Let us 
consider adiabatic processes of inserting the Aharonov- 
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FIG. 1: The Hall conductance a^^y (solid line) as a function 
of Fermi energy Ef, and density of states versus energy E 
(dashed line), at a lattice size 32 x 32. The magnetic flux 
per plaquette (</>) is 1/16, the correlation length (A) is 8, and 
disorder strength (W) is 0.01. 




FIG. 2; (Color online) Contour plot of randomness realiza- 
tion for A = 8 (a). The density probability distribution of 
several states, the localized state in the tail of lowest Landau 
band (b), the critical state in the lowest Landau band (c), 
and the critical state in the second lowest Landau band (d). 
The contour line indicated with the dashed-dotted line in (a) 
corresponds to the energy of the state of (c). The dotted line 
in (a) is the contour line with the same energy. 



Bohm (AB) flux into the holes of the torus (See Fig. [3]). 
Here it should be noted that the topology of the square 
lattice used in the present paper is a torus due to the 
periodic boundary condition. The penetration of flux $0 
(<i>i) into the hole inside (outside) of the torus, which is 
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equivalent to twisting the boundary phases 9y {Ox), could 
induce an infinitesimal electric field in the a;(?/)-direction. 
The charge-current response to the electric field is given 
by the Byers-Yang formula: 

J _ dE dE 
where E is the total energy of the system. 




FIG. 3: Two possible insertion of Aharonov-Bohm flux 
through a torus. 

When the AB flux in the unit of the magnetic flux 
quantum h/e is adiabatically turned on from to 27r, 
the spatial distribution of the probability density changes 
drastically, depending on the sign of ly as shown in Fig. 3] 
Fig. m^a) shows the energy as a function of the AB flux 
<I>o for the critical state in the lowest band. Fig. [Ifb) or 
(c) shows the typical probability density when the energy 
increases (decreases) monotonously, i.e., ly > (ly < 0). 
It turns out that this probability density corresponds to 
the percolating path along which the current flows. Ac- 
tually, Fig. mb) or (c) corresponds to the contour line 
shown as the dotted line (or the dashed-dotted line) in 
Fig. IHa). These two contour lines, which corresponds to 
the percolating paths, indicate the same energy in the 
random potential. While the $o dependence is drastic, 
we can hardly flnd a variation with $i in the total en- 
ergy and the distribution of probability density. This 
insensitivity to the flux <i>i is attributed to the fact that 
there is no percolating path in the x-direction. We also 
found the localized state is insensitive to both $o and 
$1, because there is no current path in the x direction 
nor in the y direction. The correspondence between the 
percolation path and the direction of the current can be 
understood in the view of a classic picture. In the system 
with boundaries, the electrons form skipping orbits along 
the boundaries corresponding to the cyclotron motion. 
The boundaries in this picture are formed by mountains 
of the random potential, which correspond to the dotted 
and dashed-dotted lines in Fig. [IJa). In this meaning, 
the percolation paths are thought to be the edge states. 

To take into account experimental situation, we con- 
sider averaged Hall conductance {<Jxy), which is defined 
as an average over different randomness realizations. Ef- 
fect of increasing the disorder strength W on {a^y) as a 
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FIG. 4: (a) Energy dependence on AB flux. The distribution 
of probability density of the critical state in the lowest band, 
when ly > (b) and when ly > (c). 



function of the filling factor is shown in Fig. [5ja). With 
increasing W, one can see the broadening effect on the 
steps of the Hall conductance. Moreover, the step in the 
lowest Landau band systematically shifts to the higher 
filling. This behavior is consistent with the floating pic- 
ture. On the other hand, the behavior of (axy) as a func- 
tion of energy is entirely different as shown in Fig. [5Kb). 
Sharp steps at a weak randomness W = 0.001 are broad- 
ened completely by randomness in Fig. [51(b). This is be- 
cause the randomness is not only broadening the Landau 
bands, but also shifts the position on Landau bands de- 
pending on the mean value of random potential in each 
realization. This suggests, it is difficult to detect the 
ffoat-up in energy when the sample measured in experi- 
ment has many domain walls. 

In conclusion, we have numerically calculated the Hall 
conductance, i.e., the Chern number, on a 2D system 
with a uniform magnetic field and a spatially correlated 
random potential. We have found that numerical method 
proposed by Fukui et al. is suitable to obtain the Chern 
number in a system with a random potential. We have 
also investigated the spatial distribution of the wave func- 
tions to study the topology of extended states and lo- 
calized states. The probability density of the extended 
states at the center of the Landau bands forms a cur- 
rent path associated with the percolation of equipotential 
lines. On the other hand, for the localized states in band 
tails, the charge density accumulates around the valley 
(or mountain) of the potential. The topological differ- 



4 







w =0.001 
w-0.10 
w - 0.20 




w = 0.001 

















-4.0 -3.5 -3.0 

% 



-1.5 



FIG. 5: (Color online) The averaged a^y as a function of 
filling factor (a) and Fermi energy Ef, at various strength of 
randomness (b). The correlation length (A) is 8. The assemble 
average is taken by 320 randomness realization. 



ence of wave function results in the different response to 
the AB flux: the energy of the extended state is sensitive 
to the AB flux, while the AB flux has almost no effect on 
the localized state. Finally, after taking an average of the 
Hall conductance over different randomness realizations, 
floating feature to higher filling is observed. 
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